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Abstract 

We study a class of integrable alternating {81,82) quantum spin 

chains with critical ground state properties. Our main result is the 
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description of the thermal Drude weight of the one-dimensional al- 
ternating spin chain as a function of temperature. We have identi- 
fied the thermal current of the model with alternating spins as one 
of the conserved currents underlying the integrability. This allows 
for the derivation of a finite set of non-linear integral equations for 
the thermal conductivity. Numerical solutions to the integral equa- 
tions are presented for specific cases of the spins Si and 82- In the 
low-temperature limit a universal picture evolves where the thermal 
Drude weight is proportional to temperature T and central charge c. 
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1 Introduction 



In recent years considerable progress has been achieved in the understand- 
ing of transport phenomena in low- dimensional strongly correlated quantum 
systems pQ . Thermal transport in (quasi) one- dimensional systems has been 
investigated on the theoretical as well as on the experimental side [2]. The 
existence of anomalous heat transport, indicated by for instance a non-zero 
Drude weight, has been established in particular for integrable quantum sys- 
tems. The temperature dependence of the thermal conductivities of spin-^ 
Heisenberg chain compounds was measured and revealed anomalous trans- 
port properties [21111 IS]- The experimental results point into the same direc- 
tion as the theoretical findings obtained by the Bethe ansatz technique |S] 
for the spin-| Heisenberg chain and the spin-S* Heisenberg chain [7]. 

Here we are interested in the computation of the thermal conductivi- 
ties of the general case of integrable alternating spin {81,82) chains using 
Bethe ansatz techniques. Our main result is the computation of the thermal 
Drude weight yielding a finite value for finite temperatures implying ballistic 
thermal transport. Specifically, we have obtained a finite set of non-linear 
integral equations providing the thermal conductivity as function of temper- 
ature. These equations are solved numerically for specific values of 81 and 

The paper is organized as follows. In section [21 we outline the basic 
ingredients of the Bethe ansatz techniques needed in this work. In section |3l 
we discuss the derivation of the set of non-linear integral equations for the 
description of the thermal Drude weight. Here, we also present our numerical 
and analytical findings for the solution of the non-linear integral equations. 
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Our conclusions are given in section HI 



2 Integrability and conserved currents 

In the theory of integrable models the generating function of quantum in- 
tegrals of motion is an object playing the role of a transfer matrix of the 
associated classical statistical model. The transfer matrix is the trace of an 
ordered product of Boltzmann weights defined on the square lattice. Specifi- 
cally we can define the transfer matrix T(A) of a rotational invariant classical 
vertex model associated with the alternating spin chain (^i, S2) as a product 
of two ordinary transfer matrices t^^^\\) 

T(A) = t^^'\X)t^^^\X), (1) 



t^'^X) = Tr,[£i5^^)(A)£l5!l)(A)...£g'^^)(A)£2'"^)(A)J, (2) 

where we have assumed that L is an even integer number. 

The operator C^fj'^^\X) is a rational solution to the Yang-Baxter equation 
and is obtained by for instance the fusion process [8j. It is given by 

S1+S2 

;=|Si-S2| 

where /.(A) = njm (fef) U'T^'"'"\>^ + \S2- S,\+j) and Pi is the usual 
SU{2)i projector 

P,= 'ff' (4) 

Xl — Xk 

fc=|Si-S2l 
k^l 

with Xl = \ [/(/ + 1) - + 1) - 52(^2 + 1)] and the SU{2) generators 

Sa = {S^, Sy, S^) for a = 1, 2. The operator (E]) has the following symmetry 



properties, 

Unitarity: C[l''''\X)C[l'^''\-X) = C5i,&(A)Idi (g) Ids, (5) 

Crossing: C[f'^'\\) = {-l)''^'Mijr[f'^'\-X - l)*Wf \ (6) 

Regularity: £[f'''\o) = {2S,)\Pi^'\ (7) 

where (s^s^W = U^j=^^'^'''^'\{\S2 ~ Si\ + j)^ - A^). The matrix M is an 
anti-diagonal matrix M^j — — (— 25i+2-j and Pi^2 is the permutation 
operator. Prom now on, we assume that S2 > Si without loss of generality. 
The transfer matrix T(A) can conveniently be re-written as 



T(X) = TV^ 



(8) 



L^(A)...L-r(A) 
where j — {2j — 1, 2j), a — (ai, 02) and 

= (>cS?^(A)4Sf ^(A)) ^(A)£Sf )(A)) . (9) 

The conserved currents J'^"'^ are obtained by taking logarithmic deriva- 



tives of the transfer matrix T{X) 

aA" 



^^"^ = IJ!rinT(A)[^^, (10) 



evaluated at the regular point Lgg(O) = CqP-j, where P-g = Pa^'lPa^'l and 
Co = (25i)!(2^2)! n?S('52 - Si+jf. 

The quantum Hamiltonian of the alternating spin chain is the first non- 
trivial conserved current H — J^^^\ which results in 

^ = E W' ^^^"^^ hk^i = -c/kM'^dX^kM'i^^) ,^0' 
fe=l 
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where periodic boundary conditions are assumed. For illustration, the Hamil- 
tonian for case Si = 1/2, 5*2 = 1 is given explicitly by [9] 



2 ^ / - 5 - ^ 31 ^ 1 



4 

even i ^ ^ ' 

The exact computation of the thermal conductivity and its Drude weight 
as a function of temperature is possible thanks to the fact that one can 
identify the thermal current operator with one of the conserved currents 
underlying the integrability of the model. 

In order to establish the desired connection, we consider the local con- 
servation of energy in terms of a continuity equation. This relates the time 
derivative of the local Hamiltonian h^^^ to the divergence of the thermal 
current j^, h = — Vj^. 

As the time derivative leads to the commutator with the Hamiltonian, 
we obtain 



h 



k,k+l 



where the local energy current is given by 



^k-iS' ^k,k+l 



and the total thermal current is JTe = X]fe=i • 



(13) 



(14) 



On the other hand, by inspection we find that the expression for jTs 
and the second logarithmic derivative of the transfer matrix J'^'^^ are closely 
related 



Je = J^"-^ + i 



.L52S(A) 



2 9A2 



A=0 



(15) 



where S(A) = IlKj j<2 Csi.Sjl-^)- Therefore the thermal current can be iden- 
tified as one of the conserved currents underlying the integrability. This 
property is exploited for studying the thermal Drude weight of the {81,82) 
alternating spin chain in the thermodynamical limit. 

We would like to stress that this identification is possible in the case 
of rotational invariant vertex model discussed here. In this case we have 
the staggering of spins 81 and 5*2 in the vertical and horizontal directions 
of the classical vertex model. However, a similar identification for the non- 
rotationally invariant case discussed on [7] is still an open question. 



3 Thermal Drude weight 

The transport coefficients are determined from the Kubo formula [10] in 
terms of the expectation value of the thermal current JTe, such that [21 E] 

DtniT) = (Jl) , /3:=i. (16) 

In order to calculate the expectation value we introduce an auxiliary 
partition function Z, 

Z = Tt [exp {-pn - A„ . (17) 

In this way, we obtain the expectation values of J^'^'^ through the logarithmic 
derivative of Z, 

(^)^'"^l..oHJi>. (18) 

where we used the fact that the expectation value of the thermal current in 
thermodynamical equilibrium is zero {Je) = 0. 
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To compute the partition function Z, we follow the procedure developed 
in [6] . We rewrite the partition function Z in terms of the row-to-row transfer 
matrix such that 



Z = hm Tr [T{ui) . . . T{un)T{0)-^] , 

exp I lim V{lnT(Mi) -InTfO)} 



N^oo 

Tr 



1=1 



(19) 



The numbers Ui, . . . ,un are chosen in such a way that the following relation 
is satisfied, 



N 



d 



lim V{lnT(MO -lnT(O)} = lnr(x 



dx 



x=0 

(20) 



In what follows we shall extend the results of [7] in order to include 
rotationally invariant classical vertex model associated with the alternating 
spin chain. In doing so, we need to introduce quantum transfer matrices 
rpiSi,{Si,S2}) g^j^^ rp{S2,{Si,S2}) ^Qj. gpins 5*-^ and 82- In this way, the partition 

function Z in the thermodynamical limit can be written in terms of the 
largest eigenvalues of the above mentioned quantum transfer matrices. 



lim ^InZ = -lnA(^i)(0) + - In A(^2)(0), 
L— >-oo L 2, 2 



(21) 



with 



ri— 1 



lnA(^^)(0) = (-/3 + A„,^-^)fi(a;)^^^+(ir*lny(^^))(0), (22) 
+ An|^)^2(x) 



lnA(^^)(0) 
where K{x) = 



x=0 



(is: * In 55) (0), (23) 



cosh \nx] 



and the functions Si{x) are given by Si{x) 
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/- 



dk. 



-oo 2cosh[fc/2] ' 



^ _ e-l'=l(^^+^)sinh[|A;|(5i-i)] + e-l'=l(^2-^)sinh[|A;|(^i + i)] ^^^^ 



sinh[|A;|/2] 
\k\S2 sinh[|fc|5'i] + sinh [[fclS'a] 



(25) 



sinh[|A;|/2] 

The auxiliary functions Y^^'^\x), B{x) and B{x) are required to satisfy a 
closed set of functional equations. The procedure of deriving from these 
equations a set of integral equations is similar to that described in Ref. [TT] . 
therefore we just present the final result 

/ In 2/(5) (x) \ / \ [ lnF(5)(a;) \ 



x] 



\ny^^^~^\x) 
In b{x) 
y ln6(a;) y 



+ A„£^)c/(x) 



{-(3 + Xn£^)d{x) 

V i-f3 + Xn£^)d{x) J 



lnF(5i)( 



x] 



In y (52-1) (^3,) 

ln5(x) 
In^fx) 



(26) 



where d{x) = ^^^-^^^^-^ and the symbol * denotes the convolution / * g{x) 

IZf(^-y)9iy)t 
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The kernel matrix reads explicitly 



JC{x) 



\ 





K{x) 










K{x) 




where F{x) = J 



K{x) 


K{x) 







OO Q—\k\/2+ikx 





K{x) 















K{x) 


K{x) 
K{x) 






K{x) 
F{x) 
-F(x-i) 






K{x) 

-F{x + i 

F(x, 

'(27) 



\x) J 



j^dk and K(x) was defined in Eq. 



-OO 2cosh[fc/2] 

These equations have a different driving term structure in comparison 
with those for the homogeneous case discussed in [7]. Therefore, the equa- 
tions fl2B]) constitute a new set of equations suitable for the description of 
the thermal Drude weight of the alternating spin chains. 

Finally, the thermal Drude weight can be written as 



r(2y 



_d_ 



[lnA(^^)(0) + lnA(^= 



A2=0 



(28) 



In Figure [H we show numerical results for the thermal Drude weight as 
a function of temperature for specific choices of the spins {Si,S2). At low 
temperatures the data exhibit a linear temperature dependence. The low- 
temperature asymptotics is accessible by an analytical treatment of the non- 
linear integral equations similar to [11]. From this we find that the thermal 
Drude weight Dth{T) is generally a linear function of temperature and is 
proportional to the central charge of the system c 
explicit expression reads 



gT+(f^[l2]. Its 



TiVgC 

~V1 



(29) 
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(S., ,82) case 




O.S 1 1.5 2 



Figure 1: Thermal Drude weight Dth{T) as function of temperature for 
(Si = 1/2, ^2 = 1), {Si = 1/2, S2 = 3/2), (5i = 1/2, ^2 = 2), {Si = 1, ^2 = 
3/2), (^1 = 1, ^2 = 2) and (^1 = 3/2, ^2 = 2). 



where Vg = 27r is the sound velocity. These results are in agreement with 
those of the special case of homogeneous spin chains (5*1 = ^2 = S") [7]. 

The low-temperature asymptotics (!29|) is remarkable. This result may 
be derived on the hypothetical grounds of a particle picture of the thermal 
transport. If we assume that the elementary excitations of the system are 
the carriers of the thermal transport, we would expect that the thermal con- 
ductivity ft of a system with finite mean free path I (due to imperfections of 
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the lattice) is k, = Vg c{T) I where Vg is the velocity of the elementary exci- 
tations and c(T) is the specific heat. On the other hand, the existence of a 
finite mean free path defines a mean life time r over which the Drude weight 
broadens to give a finite conductivity n = Dfh t. The condition that both 
formulas agree is Dfh = v1 c{T) which is -up to a numerical factor of order 
1- equivalent to (!29|) . It is remarkable that this reasoning yields the cor- 
rect result despite the insufficient grounds: For generic Tomonaga-Luttinger 
liquids and conformally invariant field theories there are no delta-function 
peaks in the spectral function of the single-particle Green's functions. 

4 Conclusion 

In this paper we have obtained a set of non-linear integral equations allowing 
for the explicit calculation of the thermal Drude weight for the integrable 
alternating (5*1, 6*2) spin chains. We have solved the equations numerically 
for specific values of 5*1 and 5*2. At low-temperatures we observe linear tem- 
perature dependence of the thermal Drude weight. This linear behavior is 
confirmed by the analytical low-temperature asymptotic solution showing 
that the coefficient of proportionality is the central charge of the system 
c = -|- -^^z^^^- The existence of a non-zero Drude weight at finite 
temperature is a signature of anomalous thermal transport [2] . We hope our 
findings will be useful for experimental investigations of transport properties 
of quasi one-dimensional mixed spin systems, e.g [T3]. 

Finally, we would like to remark that our results may be further gener- 
alized by including the alternation of an arbitrary number of different spins 



12 



Si,i = 1, . . . , M. Another direction of generalization would be the extension 
to anisotropic spin chains. 



Acknowledgments 

The authors thank the Volkswagen Foundation for financial support. G.A.P. 
Ribeiro also thanks FAPESP for financial support. 



References 

[1] X. Zotos, Finite temperature Drude weight of the one- dimensional spin- 
1/2 Heisenberg model, 1999 Phys. Rev. Lett. 82 1764. 

[2] X. Zotos and P. Prelovsek, Transport in one dimensional quantum 
systems (review) in "Strong interactions in low dimensions", series 
"Physics and Chemistry of Materials with Low Dimensional Struc- 
tures" , eds. D. Baeriswyl and L. Degiorgi, Kluwer Academic Publishers 



p. 347-382 (2004); [i^v:cond-mat/0304630V2 . 



[3] J. Takeya, I. Tsukada, Y. Ando, T. Masuda and K. Uchinokura, 
Thermal conductivity of Mg-doped CuGe03 at very low temperatures: 
Heat conduction by antiferromagnetic magnons, 2000 Phys. Rev. B 62 
R9260; A.V. Sologubenko, K. Gianno, H.R. Ott, A. Vietkine and A. 
Revcolevschi, Heat transport by lattice and spin excitations in the spin- 
chain compounds SrCu02 and Sr2Cu03, 2001 Phys. Rev. B 64 054412. 



13 



[4] N. Hlubek, P. Ribeiro, R. Saint-Martin, A. Revcolevschi, G. Roth, G. 
Behr, B. Bchner, C. Hess, Towards ballistic heat transport in the S=l/2 
Heisenberg chain compound SrCu02 larXiv:0908. 168II 

[5] A.V. Sologubenko, T. Lorenz, H.R. Ott, A. Freimuth, Thermal con- 
ductivity via magnetic excitations in spin-chain materials, 2007 J. Low 
Temp. Phys. 147 387. 

[6] K. Sakai and A. Kliimper, The thermal conductivity of the spin- 1/2 
XXZ chain at arbitrary temperature, 2002 J. Phys. A: Math. Gen. 35 
2173; Non-dissipative thermal transport in the massive regimes of the 
XXZ chain, 2003 J. Math. A: Math. Gen. 36 11617. 

[7] G.A.P. Ribeiro and A. Kliimper, Thermodynamics of antiferromagnetic 
alternating spin chains, 2008 Nucl. Phys. B 801 247. 

[8] P.P. Kulish, N.Y. Reshetikhin and E.K. Sklyanin, Yang-Baxter equa- 
tion and representation theory-1, 1981 Lett. Math. Phys. 5 393. 

[9] H.J. de Vega, L. Mezincescu and R.I. Nepomechie, Thermodynamics of 
integrable chains with alternating spins, 1994 Phys. Rev. B 49 13223. 

[10] R. Kubo, Statistical- Mechanical theory of irreversible process -1: Gen- 
eral theory and simple applications to magnetic and conduction prob- 
lems, 1957 J. Phys. Soc. Japan 12 570. 

[11] J. Suzuki, Spinons in magnetic chains of arbitrary spins at finite tem- 
peratures , 1999 J. Phys. A: Math. Gen. 32 2341. 



14 



[12] S.R. Aladim and M.J. Martins, Critical behaviour of integrable mixed- 
spin chains, 1993 J. Phys. A: Math. Gen. 26 L529. 

[13] O. Kahn, E. Bakalbassis, C. Mathoniere, M. Hagiwara, K. Katsumata, 
L. Ouahab, Metamagnetic Behavior of the Novel Bimetallic Ferro- 
magnetic Chain Compound MnNi(N02)4(en)2 (en=Ethylenediamine), 
1997 Inorg. Chem. 36 1530; R. Feyerherm, C. Mathoniere and O. 
Kahn Magnetic anisotropy and metamagnetic behaviour of the bimetal- 
lic chain MnNi(N02)4(en)2 (en — ethylenediamine), 2001 J. Phys.: 
Condens. Matter 13 2639. 



15 



